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Представлена методика анализа прочности метаемых сборок, которая позволяет моделировать их поведение при интенсивном нагружении 
и в некоторой мере заменять дорогостоящие экспериментальные исследования по отработке конструкции поддонов и метаемых элементов, 
подбору материалов, величин внешних нагрузок и т.д. относительно дешевым и оперативным численным экспериментом. Задача определе- 
ния упругопластического напряженно-деформированного состояния решается с помощью вариационно-разностного метода, реализующего 
вариационный принцип Лагранжа методом конечных разностей. Физические соотношения принимаются согласно теории малых упругоп- 
ластических деформаций, а геометрические соотношения берутся в виде уравнений Коши. Физически нелинейная задача решается методом 
переменных параметров упругости. Рассмотренные примеры подтверждают известную из расчетной практики возможность квазистатического 
подхода к оценке прочности подобных конструкций. 


1. Введение 

Проблемы высокоскоростного метания связаны с 
интенсивным импульсным нагружением как самих раз- 
гоняемых элементов, так и вспомогательных устройств 
- под донов. Повышение скорости метания элементов при- 
водит, как правило, к увеличению действующих на сборку 
нагрузок. Процесс метания характеризуется высокой ин- 
тенсивностью действующих на конструкцию давлений и 
перегрузок, что может привести к нежелательному уровню 
деформаций и напряжений и, как следствие, разрушению 
материала самой сборки и метательной установки. 

Экспериментальные исследования компоновочных 
схем метаемых сборок являются необходимым и обя- 
зательным этапом при их проектировании и отработке. 
Поэтому, учитывая большие временные, финансовые 
и экономические затраты при проведении этих работ, 
естественным является стремление использовать мате- 
матические методы механики деформируемого твердого 
тела (МДТТ) для моделирования поведения элементов 
конструкции и описания ее деформирования при метании 
с целью сокращения экспериментальных испытаний. 

В работе проводится анализ деформирования ме- 
таемых сборок в квазистатической постановке, хотя 
расчет напряженно-деформированного состояния (НДС) 
подобных конструкций необходимо проводить с исполь- 
зованием динамических уравнений. С точки зрения МДТТ 
процесс разгона элемента является достаточно длинным 
во времени. Но для импульсов внешней нагрузки большой 
длительности - несколько миллисекунд - численный счет 
на основе этих уравнений становится неустойчивым и 
требует неоправданных затрат времени на ЭВМ. Поэтому, 
несмотря на относительно высокую скорость протекания 
процессов при метании, в большинстве случаев можно 
использовать результаты статического расчета. В литера- 
туре, в частности, для оболочечных конструкций прини- 
мается [1] возможность использования квазистатического 
(пренебрегающего инерционными членами в уравнениях 
движения во всех или некоторых направлениях) или стати- 
ческого решения, если за время нагружения упругая волна 
успевает пробежать расстояние порядка максимального 
габарита конструкции примерно 8 ... 1 0 раз. Очевидно, еще 
в большей мере это относится к пространственным конс- 


трукциям, так как они имеют более высокие частотные ха- 
рактеристики, чем оболочечные конструкции. В конечном 
счете учет инерционных членов приводит к получению 
волновой картины распространения деформаций, которая 
становится несущественной после нескольких пробегов 
волны по конструкции, так как в результате дисперсии и 
интерференции волн постепенно устанавливается картина 
НДС, характерная для статического нагружения. 

Решение осесимметричной исходной динамической 
задачи в квазистатической постановке возможно благо- 
даря введению в расчетную схему нагруженіи массовых 
сил, являющихся следствием использованіи связанной со 
сборкой неинерциальной системы координат г, ср, г [2]. 
При метании сборка движется поступательно с ускоре- 
нием. Согласно принципу Даламбера, "... всякое тело, 
движущееся с ускорением, можно рассматривать как 
неподвижное, добавив к реально действующим силам 
фиктивные силы инерции” [3]. Возникающие осевые 
силы инерции равняются произведению ускорения 
конструкции в данный момент времени на его массу. 
Для работы с появляющимися фиктивными массовыми 
силами вводится понятие коэффициента перегрузки 
(или просто перегрузки К), который определяется отно- 
шением ускорения тела а в данный момент к ускорению 
силы тяжести §: N =а/§. Получаемая система нагрузок 
является самоуравновешенной, а из рассмотрения 
исключается постоянная составляющая осевых пере- 
мещений конструкции. Поэтому на первом этапе задачу 
можно формулировать и решать как квазистатическую. 
Возможность использования квазистатического под- 
хода для решения задачи высокоскоростного метания 
подтверждается результатами сравнения с динамичес- 
ким решением для случая умеренных пластических 
деформаций. 

Определение упругопластического НДС про- 
водится вариационно-разностным методом (ВРМ) 
[4, 5] с использованием для представления производных 
естественных аппроксимаций Ноха [6] и быстродейс- 
твующего метода Гаусса [7-9] для реализации системы 
линейных алгебраических уравнений большого порядка 
относительно искомых перемещений. Физически нели- 
нейная задача решается с помощью соотношений дефор- 
мационной теории пластичности и метода переменных 
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Рис. 1. Расчетная схема метаемой сборки 
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Рис. 2. Распределение интенсивностей деформаций е. на поверх- 
ности ОЕЭ для варианта 2 (кривая 1) и варианта 4 (кривая 
2 ) 

параметров упругости. Нелинейная зависимость интен- 
сивности напряжений от интенсивности деформаций 
аппроксимируется ломаной двухзвенной линией. 

2. Примеры расчета метаемых сборок 

Ниже рассматриваются примеры численного конс- 
труирования и анализа НДС нескольких вариантов ме- 
таемых сборок. Дабы исключить повторения, отметим 
общие моменты в постановках модельных задач для 
сборок, отличающихся друг от друга конструктивным 
решением, формой метаемого элемента и поддона, 
применяемыми материалами, величинами прилагаемых 
нагрузок и т.д. 

Задачи решаются в осесимметричной квазистатичес- 
кой постановке. Объектом исследования является тело 
вращения неканонической формы, находящееся под 
действием системы поверхностных нагрузок (давление 
газа, давление метаемого элемента) и массовых сил с ин- 



тенсивностью N. Стенки метательной установки в пер- 
вом приближении предполагается недеформируемыми, 
трением между ними и поддоном пренебрегается, хотя 
его можно учесть заданием касательных напряжений. 
Здесь же для внешней поверхности поддона задаются 
нулевые значения для радиальных перемещений. На 
оси симметрии сборки Ог касательные напряжения и 
радиальные перемещения равны нулю. На свободных 
поверхностях полагаются равными нулю нормальные и 
касательные напряжения. На контактных поверхностях 
конструктивных элементов сложных поддонов и мета- 
емых элементов задаются условия равенства перемеще- 
ний и напряжений, т.е. условия непрерывного контакта. 
Вследствие симметрии конструкций относительно оси 
вращения г на рисунках представлены лишь половины 
осевых сечений рассматриваемых сборок. 

Одна из возможных схем метаемой сборки приведе- 
на на рис. 1 : это поддон и метаемый элемент - шарик. 
Шарик считается недеформируемым. Поддон конс- 
труктивно состоит из двух частей - дюралюминиевого 
цилиндра (заштрихованная часть), в котором наполо- 
вину утоплен шарик, и непосредственно пластикового 
поддона. Устройство разгоняется газом с давлением Р р 
приложенным к левому торцу поддона ОА. На поверх- 
ности БЕ цилиндра приложено в осевом направлении 
давление Р ш от действия метаемого шарика. На сборку 
действуют также и массовые силы с интенсивностью N. 
направленные вдоль оси Ог противоположно давлению 
Р г . Геометрические размеры поддона и шарика брались 
следующие: длина поддона АВ = 0,8 см, ОН = 0,5 см, 
К.! = 0,4 см, К 2 = 0,3 см. К, = 0,15 см. Масса мета- 
емого шарика в зависимости от его материала - сталь 
или сплав ВНЖ - 0,1 и 0,27 г соответственно. 

Были рассмотрены четыре варианта этой мо- 
дельной задачи, которые различались величинами 
нагрузок и материалом метаемого шарика (см. 
рис. 3). Варианты 1 и 2 посчитаны для шарика, мате- 
риалом которого является сталь, а варианты 3 и 4 - для 
сплава ВНЖ. Распределение интенсивности деформа- 
ций е. на поверхности 0ЕБ для нагружения по второму 
(кривая 1) и четвертому варианту (кривая 2) представле- 
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Рис. 3. Зоны пластических деформаций сборки для разных значений давления газа Р и плотности р материала метаемого шарика 
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ны на рис. 2. Увеличение плотности материала шарика 
приводит к значительному росту е г 

На рис. 3, а изображены зоны пластических дефор- 
маций (заштрихованная область) пластикового поддона 
и дюралюминиевого цилиндра для первого варианта 
нагружения. Наибольшие значения интенсивности 
деформации е. и напряжений а. имеют место на оси 
вращения Ог под метаемым шариком, на контактной с 
шариком поверхности ЕБ, а также в цилиндре на по- 
верхности МН. Нижняя часть поддона близка к выходу 
в пластическую область деформирования. Наибольшее 
значение для цилиндра интенсивности е ішах = 3,6 %. 
Увеличение давления газа до 600 МПа (вариант 2) 
приводит к увеличению размеров зоны пластической 
деформации в поддоне (рис. 3, б) и величины интенсив- 
ности деформаций - е ітах * 6 %, а также всех остальных 
параметров НДС. 

На рис. 3, в, г (варианты 3, 4) представлены резуль- 
таты для случая, когда материал шарика - сплав ВНЖ 
с большой плотностью. Использование этого материала 
приводит к увеличению значений параметров НДС 
элементов поддона и зон пластических деформаций. 
Наибольшие значения е. для поддона и цилиндра состав- 
ляют 10, 16 и 13, 20 % для вариантов 3 , 4 соответственно . 
Но величина интенсивности массовых сил N (пере- 
грузка) с увеличением массы шарика при неизменной 
величине давления газа уменьшается, и соответственно 
уменьшается скорость метания. В описанном случае 
уменьшение перегрузки составляет а 18 %. Наиболее 
напряженным местом в поддоне является участок на 
контактной поверхности ЕБ с координатой г а 0,7 см. 

Следующей была рассмотрена задача определения 
НДС части поддона, расчетная схема которого пред- 
ставлена на рис. 4. Поддон разгоняется газом (Р = 
350 МПа), действующим на поверхность дна ОБ. На 
поверхности СО приложены давления Р от метаемого 
стержня и Р т от другого поддона, который не рассмат- 
ривается. К поддону приложены также и массовые 
силы с перегрузкой К, направленной вдоль оси 0г про- 
тивоположно давлению газа Р г Конструктивно поддон 


состоит из стального подпятника 0БЕРНС, на который 
опирается вольфрамовый стержень, и помещенного на 
подпятнике текстолитового цилиндра ЕОНР (на рис. 
4 заштрихованная область). Геометрические размеры 
брались следующими: К.^1,165 см, К,= 0,4 см, длина 
ОС = =2,0 см. 

Был проведен численный эксперимент по выяв- 
лению наиболее выгодного с точки зрения прочности 
варианта конструкции поддона путем задания положе- 
ния контактной поверхности ЕРН при фиксированной 
длине ОС. Результаты этих расчетов для пяти вариантов 
конструкции поддона представленны в таблице и на 
рис. 5, где штриховкой изображены зоны пластической 
деформации материала подпятника и цилиндра. 

Данные таблицы и рис. 5 позволяют выделить вари- 
анты 1 и 4, в которых меньше величины интенсивностей 
е. и а., а также размеры зон пластического деформи- 
рования подпятника, как наименее напряженные. Это 
объясняется более мощной конструкцией подпятника, 
нежели в вариантах 2, 3, 5. Если в варианте 1 увели- 
чить длину поддона ОС до 2,5 см, не изменяя размера 
текстолитового цилиндра, то максимальные значения 
е ; и а уменьшаются на 10... 15 %. При этом уменьша- 
ется также и величина осевой перегрузки - N=-7324, 
т.е. и скорость метания. Таким образом, для заданных 
величин давления газа и рассмотренных материалов 
для стержня, подпятника и цилиндра варьированием 
геометрии поддона снизить уровень деформаций не 
представляется возможным. 

На рис. 6 представлены схема метаемой сборки 

(а) и изменение нагрузки в зависимости от времени 

(б) . Сборка состоит из дюралюминиевого поддона 
типа “стакан” и стального элемента в виде цилин- 
дрического стержня. Геометрические параметры 
следующие: толщина днища поддона Ц ] = 2,0 см, 
длина поддона Б, = 10,0 см, длина конструкции 
Б 3 = 12,0 см, радиус стержня К 0 = 0,7 см, радиус под- 
дона К. = 1,2 см. Характеристики формы импульса 
равны т 1 = 50 мкс, т 2 = 300 мкс, т, = 350 мкс. Дефор- 
мация начала текучести материала поддона равна 
е = 0,005. 

При рассмотрении этой задачи реализована возмож- 
ность двух подходов - квазистатическо го и динамичес- 
кого, что позволяет провести сравнение полученных 
решений и оценить возможности каждого из них по 
адекватному отображению протекающих в конструктив- 
ных элементах сборки процессам деформирования. 

С точки зрения динамического подхода данный 
диапазон времени нагружения является “неудобным” 
в вычислительном плане, так как приводит к довольно 
большим затратам машинного времени при счете про- 


Таблица. Результаты численного эксперимента 
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Рис. 5. Зоны пластических деформаций для пяти вариантов конс- 
трукции поддона 
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Рис. 6. Расчетная схема метаемой сборки (а) и изменение нагрузки 
в зависимости от времени (б) 
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Рис. 7. Зоны пластических деформаций в сборке для трех вари- 
антов нагружения в случае квазистатического приложения 
давления газа 


цесса деформирования по явным конечно-разностным 
схемам. Применение этого подхода может быть оправ- 
дано в случае невозможности пренебречь волновыми 
процессами в сборке, при наличии больших упругоп- 
ластических течений и при сложном нагружении. Здесь 
используются соотношения, описывающие процесс 
упруго пластического деформирования тела вращения и 
базирующиеся на законах сохранения массы, импульса 
и энергии. Поведение материала подчиняется соотноше- 
ниям пластического течения типа Прандтля-Рейсса [10]. 
Конечно-разностная схема второго порядка аппроксима- 
ции записывается для лагранжевой сетки с использова- 
нием интегрального представления пространственных 
производных [6], что позволяет использовать нерегуляр- 
ные и деформированные разностные сетки. Для расчета 
упругопластических течений используется процедура 
приведения напряжений на круг текучести [10]. 

На рис. 7 приведены результаты расчета зон плас- 
тических деформаций в случае квазистатического 
приложения давления Р и соответствующих им мас- 
совых сил с интенсивностью № а) 693 МПа, -15000, 
б) 808,4 МПа, -17500, в) 924 МПа, -20000. Деформация 
начала текучести материала стержня равна е 5 = 0,005. 
На рисунках видно, что увеличение нагрузки приводит к 
формированию и росту в поддоне зон пластических де- 
формаций. По мере увеличения давления они образуют- 
ся вначале под торцом стержня, затем распространяются 
на всю толщину днища, и далее пластическим дефор- 
мированием охвачена почти вся нижняя часть поддона. 
Материал стержня работает в упругой области. 

Представленные на рис. 8 результаты получены для 
нагрузок Р =924 МПа, N = -20000 и меньшего значения 
величины предела текучести е 5 = 0,003 материала 
стержня, что приводит к появлению в нем пластичес- 
ких деформаций. Сравнение результатов расчета по 
квазистатической теории (рис. 8, а) с результатами 
расчета по динамической схеме при условии непрерыв- 
ного контакта между стержнем и поддоном на боковой 
поверхности (рис. 8, б) и при условиях скольжения на 
этой поверхности (рис. 8, в) приводит к выводу о том, 
что качественно картина пластического деформирования 
хорошо совпадает для первых двух случаев. 

Следует отметить, что во всех просчитанных задачах 
все напряжения - сжимающие и, следовательно, не мо- 
гут представлять опасности с точки зрения целостности 
конструкции. Деформации же имеют разные знаки, 
что и приводит к большим значениям интенсивности 
деформаций в;, величина которой определяется главным 
образом сжимающей осевой деформацией. 

Вообще говоря, эти модельные задачи могут ста- 
виться шире, нежели только как задачи определения 
НДС поддона. Необходимо получить такую конструк- 
цию, чтобы материалы конструктивных элементов 
поддона выдержали нагрузки и не разрушились при 
метании. Полученные результаты дают представление 
о характере распределения напряжений и деформаций в 
элементах поддонов для достаточно большого диапазона 
изменения давления газа и интенсивности массовых сил 
(скоростей метания), различных материалов и вариантов 
конструкций как поддонов, так и метаемых элементов, 
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Рис. 8. Зоны пластических деформаций в сборке для 
Р=924 МПа, N =-20000 в случае: квазистатического (а) и 
динамического (б) приложения давления газа при условии 
непрерывного контакта между стержнем и поддоном на 
боковой поверхности; для динамического нагружения при 
условии скольжения на боковой поверхности (в) 


что позволяет вносить необходимые конструктивные 
изменения в элементы сборки с целью уменьшения 
значений тех или иных параметров НДС. При создании 
подобных конструкций, очевидно, следует пойти путем 
использования материалов с более высокими значе- 
ниями величин физико-механических характеристик. 
Эти материалы должны иметь невысокие плотности и 
выдерживать большие деформации не разрушаясь. Од- 
новременно следует проводить исследование влияния 
геометрии конструктивных элементов поддонов с целью 
уменьшения величины параметров НДС сборки. 

Автор приносит благодарность А.В. Герасимову за 
предоставленные результаты динамического расчета 
метаемого элемента с поддоном типа “стакан”. 

3. Выводы 

По результатам анализа НДС метаемых сборок 


можно сделать некоторые выводы. Созданы квазиста- 
тическая физическая и математическая модели и реали- 
зующие их пакеты прикладных программ для расчета 
деформирования неоднородных метаемых сборок при 
нагружении их комбинацией поверхностных и массовых 
сил большой интенсивности. Известная из расчетной 
практики возможность квазистатического подхода к 
оценке прочности большинства подобных конструкций 
подтверждается в рассмотренных примерах. Это связано 
с согласованием характерного времени нагружения и 
периода пробега упругих возмущений по конструкции 
малых габаритов. Уменьшение времени расчета при 
квазистатическом подходе ( 1 мин) по сравнению с фи- 
зически более оправданным волновым подходом (120 
мин), сохраняя основные черты волнового процесса 
деформирования, позволяет выйти на решение задач 
оптимизации с целью уменьшения массы поддона без 
потери его прочности. 
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